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Sciences 



Under earth surface conditions, in ocean and natural water, calcium carbonate is ubiquitous, forming 
anhydrous and hydrous minerals. These hydrous phases are of considerable interest for their role as 
precursors to stable carbonate minerals. Atomistic simulation techniques have been employed here to 
perform a comprehensive and quantitative study of the structural and energetic stability of dry and hydrous 
surfaces of calcium carbonate polymorphs using two recently developed forcefields. Results show that the 
dry forms are prone to ductility; while hydrous phases are found to be brittle. The (001) surface of 
monohydrocalcite appears to be the most stable (0.99 J/m 2 ) whereas for the ikaite phase, the (001) surface is 
the most stable. The corresponding value is 0.2 J/m 2 , i.e. even lower than the surface energy of the 
(1014) surface of calcite. Beautiful computed morphology pictures are obtained with Xiao's model and are 
very similar to the observed SEM images. 

Carbonates are among the most abundant minerals in the environment. We can found them in aquatic 
systems forming huge amounts of scale and ocean sediment 1 , in the atmosphere 2 and as part of the 
sediment and rock record 3 . These minerals affect the chemistry of aquatic systems by regulating pH 
and alkalinity through dissolution/precipitation equilibrium. Moreover, they can govern the mobility and cycling 
of metal contaminants and radionuclides via ion exchange, adsorption, and co -precipitation reactions 4 . Calcium 
carbonate (CaC0 3 ) is one of the constituents of argillite used for the wastes disposal. It has an important role in 
restraining the swelling properties of smectitic layers and produce hydraulic and mechanical reinforcement of the 
clay phase. It is one of the most abundant materials found in sedimentary rock in all parts of the Earth surface (it 
makes up 4% of the earth's crust). CaC0 3 occurs naturally in six different forms: three crystalline polymorphs, 
calcite, aragonite, and the metastable vaterite; two hydrate phases (monohydrocalcite (CaC0 3 -H 2 0) and Ikaite 
(CaC0 3 -6H 2 0)); and amorphous calcium carbonate (ACC). The monohydrocalcite structure has 8-fold coor- 
dinated calcium ions. Water is situated between the screw axes (parallel to the c-axis), bound to the carbonate ions 
via hydrogen bridges 5 . Additionally, monohydrocalcite is found to be metastable under all conditions 6 . However, 
it does form in some predominantly freshwater environments 7 . The structure of ikaite is characterized by 
hydrated calcium carbonate chains, aligned along the a- axis and are linked together by hydrogen bonds; the 
calcium ions are 8 -fold coordinated 8 . It is found that ikaite is stable at high pressure and low temperatures in an 
alkaline aqueous environment 9 " 11 . It can be found in bottom waters in fjords and on deep sea floors 12 . 

Calcite is one of the most stable phase at atmospheric pressure and temperature, while aragonite is found to be 
the most stable polymorph in high-pressure environments, although it is possible for aragonite to form at 
atmospheric pressure in certain conditions, such as hot springs 13 . Vaterite, the metastable form of calcium 
carbonate, changes to form another polymorph, or is stabilized by water adsorption on its surfaces 14 . We can 
explain the differences between the three polymorphs from the different arrangement of the C0 3 ions in their 
crystal structures. 

Many forcefields have been developed to describe the CaC0 3 interactions 14 " 22 . Other empirical models have 
been used to explore the bulk and surface properties of CaC0 3 at ambient conditions 18 23 25 . Parker et al. 26 , using 
an atomistic model, examined the surface precipitation and dissolution processes. On the other hand, Fisler et al. 22 
used the shell model to simulate the rhombohedral carbonates and their point defects. Archer et al. 17 developed a 
new empirical potential to study the first-order phase transition between calcite I and calcite II. Recently, 
new models have been developed describing the interaction of the system with water, and even organics, with 
a view to addressing biomineralization. Raiteri et al. 20 proposed a new potential model, which gives successful 



SCIENTIFIC REPORTS | 3 : 1587 | DOI: 1 0.1 038/srep0 1587 



1 



atomic-resolution simulations of the growth of amorphous calcium 
carbonate. However, the forcefield reproduces well structural and 
thermodynamic properties of CaC0 3 but not the elastic behavior. 
The latest forcefield published by Xiao et al. 21 for CaC0 3 derived a 
new set of forcefield parameters, which reproduces accurately the 
elastic and shear modulus for bulk aragonite crystal. 

Results 

The broad objective in this paper is to use, for the first time, two latest 
forcefields developed recently by Raiteri 20 and Xiao 21 in order to find 
out the best model to analyze the surface of dry and hydrous calcite 
and its polymorphs. X-Rays diffraction patterns will be evaluated and 
compared to experimental findings. Structural and mechanical 
quantities will be evaluated for calcite, aragonite, vaterite, monohy- 
drocalcite and Ikaite phases. We will discuss about the structural and 
energetic stability of CaC0 3 polymorphs through the calculation of 
the surface and attachment energies as well as the equilibrium 
morphologies. 

Bulk properties of calcium carbonate polymorphs. Calculations 
are performed using the General Utility Lattice Program 4.0 
(GULP) 27 . Internal energy minimizations have been done in NPT 
ensemble allowing all individual ionic coordinates and lattice para- 
meters to vary. The unit cell was fully optimized in order to obtain the 
equilibrium crystal structure. The search of the local minima adopts 
the Newton -Raphson procedure, with the Broyden-Fletcher-Gold- 
farb-Shanno (BFGS) scheme to update the Hessian. Structural 



parameters of the known anhydrous (calcite, aragonite, vaterite) as 
well as hydrous polymorphs of CaC0 3 (monohydrocalcite and Ikai- 
te), have been evaluated with both models (Raiteri' s and Xiao's). 
They are compared to experimental data taken from the American 
Mineralogist database. The absolute errors of the lattice parameters a, 
b, c, and equilibrium volume are displayed in Figure 1 (a, b, c, d). We 
can see that for the case of calcite, the Raiteri's model reproduces 
better the structure with a lower error (less than 2%). However, when 
studying aragonite or vaterite structures and the hydrous phases 
(monohydrocalcite and Ikaite), our results show that the Xiao's 
forcefield reproduces well the structure of hydrated CaC0 3 mine- 
rals, since the deviations of the unit cell parameters are less than 3%. 
Figure 1 (d) illustrates that the volume evaluated with Xiao model is 
obtained with an error less than 4.5% for the dry and hydrated 
phases, while Raiteri model gives a higher error (7-8%) for mono- 
hydrocalcite and ikaite and around 4% for the dry phases. 

In Figure 2, we plot our calculated X-Ray diffraction patterns 
(k = 1.54A for calcite carbonate polymorphs using Raiteri's and 
Xiao's forcefield, compared with experimental findings. As one 
would expect, the Raiteri's forcefield reproduces perfectly the calcite 
phase than Xiao's model (Figure 2 (a)). This confirms our structural 
data as discussed above. We can see that from 20 = 70°, the curve 
obtained from Xiao forcefield does not match the experimental find- 
ings. For the case of aragonite phase (Figure 2 (b)), the obtained X- 
Rays patterns from Xiao model are very similar to the experimental 
curve. We note a shift of the Raiteri's curve along x, which is more 
pronounced from 20 = 30° to 42°. For the case of vaterite, the 
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Figure 2 | X-Rays diffraction patterns for dry and hydrous CaC0 3 polymorphs calculated with Raiteri's and Xiao's forcefields and compared to 
experimental findings: (a) calcite, (b) aragonite, (c) vaterite, (d) monohydrocalcite, (e) Ikaite. 
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Figure 3 | Energetic stability of CaC0 3 polymorphs. 



monohydrocalcite and ikaite phases, X-Rays diffraction patterns 
obtained from Xiao's model are similar to experimental curves 
(Figure 2 (c), (d), (e)). In Figure 3, the energetic stability of CaC0 3 
polymorphs is presented compared to experimental measurements 28 . 
We can see that ikaite has the lowest total energy ( — 33.9393 eV). 
Calcite (-29.9371 eV) and aragonite (-29.9358 eV) phases have a 
small energy difference, about 0.0013 eV. Vaterite with a total energy 
equal to —29.8682 eV is the less stable phase. 

Concerning the mechanical properties, the elastic constants tensor 
is calculated from the second derivative matrix at the optimized 
structures and is listed in Table 1(a) for the hydrated phases. In a 
recent work of Xiao et al. the elastic data, for the dry phases, were 
reported and compared with experiments 21 ' 29 ' 30 . We can notice that 
Cn for Ikaite is smaller than C 2 2 and C 33 what means that this 
structure is more compressive in the direction perpendicular to the 
layers (layers parallel to (100)). We have also calculated the isotropic 
bulk (B) and shear (G) moduli obtained from the Voigt-Reuss-Hill 
averages 31 . We can define B as the amount of pressure required for a 
unit volume change of the material; G measures the shear forces 
required to induce a unit shape change measured, while E is the force 
per unit of cross-sectional area needed for a unit dimensional change. 



Table la | Elastic data for monohydrocalcite and Ikaite 

Monohydro calcite Ikaite 
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C,,(GPa) 


99.3 


98.8 


41.7 


59.2 


C 22 (GPa) 






96.1 


75.4 


C 33 (GPa) 


91.9 


88.1 


82.0 


62.7 


C 12 (GPa) 


25.6 


28.9 


18.8 


20.4 


C 13 (GPa) 


34.8 


42.2 


27.8 


26.3 


C 23 (GPa) 






18.7 


16.8 


C 44 (GPa) 


52.8 


32.1 


8.0 


4.2 


C 55 (GPa) 






8.9 


2.7 


C 66 (GPa) 


36.8 


34.9 


16.2 


12.9 


a Our calculations using 


Raiteri's model, 


3 using Xiao's model. 







In Table 1(b) , we have listed the bulk modulus, B, shear modulus, G, 
Young's modulus, E and B/G for all phases. Using the ratio B/G (the 
resistance to fracture (B) relative to the plastic deformation resistance 
(G)) introduced by Pugh 32 , we can analyze the ductility of a material. 
A value of B/G higher than 1.75 indicates a tendency for ductility, 
while a value of B/G lower than 1.75 indicates a tendency for brittle- 
ness. The dry forms of calcium carbonate yield values higher than 
1.75, suggesting that calcium carbonate materials are prone to duct- 
ility. Moreover, the value of B/G on going from CaC0 3 to CaC0 3 - 
H 2 0 and CaC0 3 -6H 2 0 decreases, which suggests that the inclusion 
of water in the crystal makes the material more brittle. In fact, the 
bulk modulus decreases with about 36% from calcite to monohydro- 
calcite and about 56-82% from calcite to ikaite. 

Surface analysis. Starting from the experimental structure (a = b = 
4.988 A, c = 17.061 A; a = ft = 90° and y = 120° the optimized bulk 
structure of calcite was cut to obtain the surfaces (10l4), (1120), 
(10T0), (10Tl}Ca +2 , (lOll) C0 3 - 2 , (0U2) (see Table 2). Cleavage 
along the {10 14} plane gives rise to the characteristic rhombohedra 
shape of calcite crystals. This plane contains both Ca +2 and C0 3 " 2 
ions, making the surface charge neutral (Figure 4 (a)), and thus stable 
with a low surface energy. Furthermore, it has also a higher density of 
ions compared to other possible neutral planes. For instance, the 
(0112) and the (1011) planes are polar as they are terminated by 
either a layer of calcium or carbonate ions, leading to a positively 
or negatively charged surface respectively (Figure 4 (b), (c)). We plot 
in the top of the slabs the iso-surfaces. These are obtained using the 
Gaussian description of molecular shape. It is similar to the idea of 
rolling a probe atom over the surface of the system being studied. 

For the aragonite phase, the crystal was cut to obtain the following 
low index surfaces: (010), (100), (001), (110), (011), (101), and (111). 
The (010) and (110) surfaces are cleavage planes while the (110) is 
also a twinning plane (see Table 3). The surfaces can be cut in more 
than one way, giving rise to planes terminated by either calcium ions 
or oxygen atoms of the carbonate groups. In the case of the meta- 
stable vaterite phase, all surfaces have both calcium and carbonate 
terminated planes as presented in Table 4. Surfaces are allowed to 
relax, the calculated surface and attachment energies with both 
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Table 1 b | Elastic data for monohydrocalcite and Ikaite 

Calcite Aragonite Vaterite Monohydro calcite Ikaite 
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B 


83.1 


85.7 


80.6 


72.6 


67.7 


67.1 


88.1 


69.1 


53.4 


44.4 


36.3 


15.6 


G 


32.6 


35.7 


32.7 


38.8 


35.8 


35.8 


36.9 


18.8 


36.9 


27.6 


31.1 


13.7 


B/G 


2.5 


2.4 


2.4 


1.8 


1.9 


1.8 


2.3 


3.6 


1.4 


1.6 


1.16 


1.13 


E x 


108.6 


78.1 




185.7 


132.7 




-75.0 


58.3 


80.4 


75.1 


38.4 


45.4 


Ey 


108.6 


78.1 




81.0 


69.3 




-75.0 


58.33 


80.4 


75.1 


96.4 


66.8 


E z 


54.62 


39.45 




58.56 


76.95 




145.9 


110.6 


72.52 


60.2 


81.9 


44.2 



a Our calculations using Raiteri's model, b using Xiao's model; c33 , d : 
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Figure 4 | (a) (1014) surface of calcite; on the top are presented iso-surfaces. (b) (0112) surface of calcite with Ca termination, (c) (1011) surface of calcite 
with C0 3 termination. Blue atoms represent calcium, green atoms are carbon and red atoms represent oxygen. 



Table 2 Relaxed surface and attachment energy of Calcite 






Surface energy (J. m 2 ) 






Attachment energy (eV) 




Surface (Hexagonal 














indices) 


Raiteri's forcef 


eld Xiao's forcefield 


Other results a 


Raiteri's forcefield Xiao's forcefield 


Other results a 


(10T4) 


0.71 


0.51 


0.59 


-11.34 


-8.69 


-7.82 


(1120) 


1.1 1 


0.76 


1.39 


-33.91 


-24.58 


-30.24 


(10T0) 


1.04 


0.64 


0.97 


-33.29 


-14.37 


-78.86 


(10Tl)Ca 


1.61 


1.63 


1.23 


-11.92 


-30.63 


-31.89 


(10Tl)CO 3 


1.36 


0.97 


1.14 


-17.88 


-24.58 


-28.72 


(01T2) 


4.85 


1.56 




-206.70 


-100.78 




a 14 
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Table 3 | Relaxed surface and attachment energy of Aragonite 



Surface 




Surface energy (J.m 2 ) 






Attachment energy (eV) 




Raiteri's forcefield 


Xiao's forcefield 


Other results 0 


Raiteri's forcefield Xiao's forcefield 


Other results a 


(01 0\ Cn 


1.14 


0 RO 


0 cm 

U.7U 


-4.62 


-3.20 


- 1 0.91 


(010) COo 


1 .57 


1.10 


1 .50 


-11.96 


-10.48 


-32.68 


(1 00) 


1 .71 


1 .23 


1 .50 


-15.95 


-10.75 


-39.58 


(001 ) Cn 


1 .59 


1 .36 


1 .05 


-12.73 


-9.58 


-35.41 


(ooi \ rOo 


1 01 


0 7fi 

\J./ o 


0 as 


-5.35 


-3.92 


- 1 2.46 


n i ni Cn 


1 fH 


0 7^ 


0 


-8.61 


-6.10 


— 1 0 ^50 


(110) co 3 


1.30 


0.92 


1.04 


-28.14 


4.37 


-31.45 


(101)Ca 


1.19 


0.90 


0.99 


-10.54 


-6.18 


-15.1 1 


(101)CO 3 


1.31 


1.00 


1.08 


-14.52 


-12.03 


-16.14 


(01 1)C0 3 


0.90 


0.84 


0.69 


-5.97 


-9.08 


-6.38 


(01 1)C0 3 


1.09 


1.33 


1.16 


-14.33 


-16.37 


- 1 7.40 


(011)CO 3 


1.33 


1.20 


0.99 


-10.08 


-16.23 


-10.61 


(1 1 1) Ca 


1.30 


0.94 


1.40 


-23.39 


-3.82 


-33.19 


(in)co 3 


1.09 


0.83 


0.84 


-14.94 


-6.50 


-18.38 



Table 4 Relaxed 


surface and attachment energy of vaterite 










Surface 




Surface energy (J.m 2 ) 




Attachment energy (eV) 




Raiteri's forcefield 


Xiao's forcefield 


Other studies a 


Raiteri's forcefield 


Xiao's forcefield 


Other studies a 


(010) Ca 


1.51 


1.34 


1.35 


-92.65 


4.04 


-35.58 


(010) co 3 


0.75 


0.75 


0.62 


-37.48 


-37.48 


-8.29 


(100) Ca 


1.55 


1.34 


1.31 


-1 15.51 


4.04 


-17.50 


(100) co 3 


0.75 


1.54 


1.39 


-37.48 


-91.41 


-14.23 


(110) Ca 


1.54 


0.87 


1.18 


-146.66 


-107.35 


-25.94 


(110) co 3 


0.87 


0.98 


1.04 


-49.33 


-55.29 


-15.47 


(Oll)Ca 


1.29 


1.16 


1.22 


-41.10 


-32.70 


-59.90 


(011) co 3 


1.12 


1.12 


0.93 


-55.23 


-55.23 


-19.37 


(101)Ca 


1.35 


0.86 


1.47 


-34.66 


-20.14 


-19.33 


(101)CO 3 


1.22 


0.79 


1.07 


-43.46 


-16.62 


-18.29 


(1 1 l)Ca 


1.48 


0.78 


0.82 


-190.73 


-109.70 


-19.47 


(in)co 3 


2.39 


0.86 


0.85 


-125.23 


-100.57 


-25.47 





forcefields are given in Tables 2, 3, and 4 and are compared to other 
theoretical results 14 . 

Discussion 

In this part we will discuss about the best forcefield for modeling the 
surface of calcite and its polymorphs. We will analyze, for the first 
time the equilibrium morphology of dry and hydrous phases of 
CaC0 3 . 

In the case of calcite (Table 2), it is clear that the relaxed (1014) 
surface is the most stable surface, with a lowest surface value and 
attachment energy. The corresponding values are found to be equal 
to 0.51 J/m 2 and — 8.69 eV. We notice also that Xiao's forcefiled gives 
lower values, which agree better with the experiment and other 



Table 5 | Relaxed surface and attachment energy of monohydro- 
calcite 

Surface energy (J.m" 2 ) Attachment energy (eV) 

(hkl) Raiteri's model Xiao's model Raiteri's model Xiao's model 

(100) 1.48 1.54 -13.75 -7.82 

(010) 1.39 1.21 -10.94 -6.62 

(001) 1.00 0.99 -24.94 -5.66 



theoretical data 33 " 35 . The other polar surfaces are unstable by them- 
selves and will either reconstruct to neutralize the dipole moment or 
adsorb ions to remove the dipole moment. However, they can be 
stabilized when grown on charged substrates. 

Concerning the surface of the aragonite phase and according to the 
above structural calculations, we remark (Table 3) that the surface 
energies obtained with Xiao's model agree better with the other 
theoretical works 14 . In the other hand, we can see that the surface 
energy of the pure (010) surface quoted in Table 3 shows the calcium 
terminated plane to have a lower surface energy than the carbonate 
surface and hence to be the most stable surface plane. This is due to 
the rotation of carbonate groups rotate, which tend to lie flat in the 
surface rather than upright as in a bulk termination. On the other 
hand the calcium-terminated plane shows hardly any relaxation. In 



Table 6 | Relaxed surface and attachment energy of Ikaite 

Surface energy p.m" 2 ) Attachment energy (eV) 

(hkl) Raiteri's model Xiao's model Raiteri's model Xiao's model 

(100) 1.46 0.52 -13.65 -8.19 

(010) 0.35 0.37 -7.48 -5.79 

(001) 0.23 0.20 -3.04 -2.57 
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addition, the attachment energies agree well with the surface ener- 
gies, where we notice that the calcium plane is favored over the 
carbonate plane through its low value indicating a small growth rate. 
The same remarks for the (110) surface. Concerning the (Oil) 
surface, it has different cuts with one or more surface oxygen ions. 
However, their surface energies vary greatly as shown in Table 3 from 
a very stable surface with y = 0.84 J m" 2 to the least stable with a 
surface energy of 1.33 J m~ 2 . We notice that the attachment energies 
show the same trend as the relaxed surface energies with a very low 
value for the most stable plane. Consequently, this surface would be 
expressed in the growth morphology. 

For the relaxed surface of vaterite phase, we note that our results 
based on Xiao's model are generally in good agreement with the 
theoretical calculations 14 . It is clear from Table 4 that the (010) car- 
bonate plane is the dominant surface, with the lowest surface and 
attachment energies of the series. Moreover, the (100), (001), and 
(110) surfaces have relatively large surface and attachment energies, 
even if some surfaces have been relaxed sufficiently. We remark that 
the (011) calcium plane and both (101) surfaces are dipolar and both 
need extensive relaxation on energy minimization. 




Figure 6 | Equilibrium morphology of aragonite compared to SEM ii 



In Table 5 and 6, we present our results for the relaxed surface 
of hydrated CaC0 3 , monohydrocalcite and Ikaite, using Raiteri's 
and Xiao's forcefields. For Monohydrocalcite (CaC0 3 -H 2 0), the 
surface energies found using Xiao' model are 1.54, 1.21 and 0.99 J/ 
m 2 for the x, y, and z surfaces, respectively (Table 5). The corres- 
ponding attachment energies are —7.82, —6.62 and —5.66 eV for 
the x, y, and z surfaces, respectively. These results show that the 
(100) surface is less stable than the other two surfaces, which 
implies that the (100) surface will not appear in the monohydro- 
calcite equilibrium morphology. For the case of relaxed Ikaite 
surface, the obtained surface energies (Table 6) are equal to 
0.52, 0.37, 0.20 J/m 2 for (100), (010) and (001) planes, respect- 
ively. The corresponding attachment energies are —8.19, —5.79, 
-2.57 eV for (100), (010) and (001) planes, respectively. We can 
clearly see that the (001) surface is the most stable (0.2 J/m 2 ), with 
a value lower than that found for the most stable (1014) 
surface of calcite However, the (100) plane with 0.52 J/m 2 would 
not be expressed in the growth morphology. 

In order to measure the relative stability of the surfaces, it is neces- 
sary to compute the equilibrium morphology of a crystal, which is 
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Figure 7 | Equilibrium morphology of vaterite compared to SEM images : 



determined by the surface energy and the related growth rate for the 
various surfaces. The Gibbs-Wulff theory provides a simple math- 
ematical approach that predicts the equilibrium crystal shape, which 
should possess a minimal total surface energy for a given volume 36 ' 37 . 
This induces that a surface with a high surface free energy has a large 
growth rate, and this fast growing surface will not be expressed in the 
equilibrium morphology of the resulting crystal. Therefore, we con- 
sider only surfaces with low surface energies and hence slow growth 
will be expressed. It is well known that at 0 K, the surface free energy 
is a close approximation of the surface energy as calculated by static 
lattice simulations. In fact, this is due to the entropy term, which is 
found small when included in the surface free energy. Thus, the 
surface energies can be assumed to determine the equilibrium mor- 
phology of the crystal. 

Figures 5, 6 and 7 show the equilibrium morphologies of dry 
calcium carbonate (calcite, aragonite, vaterite) obtained for all 
allowed faces and compared to the experimental images of scanning 
electron microscopy (SEM) 38 . We remark that our computed 
morphologies are very similar to the experimental ones. For the 
hydrated phases, we show in Figure 8 (a) the equilibrium morpho- 
logy of Ikaite compared to the SEM image 39 . In Figure 8(b) the pre- 
dicted equilibrium morphology for monohydrocalcite phase is 
presented since no experimental images are available. These results 
confirm once again the potential models employed and support 
strongly our findings. 




Methods 

Potential model and surface simulation. In Raiteri's model 19 , the oxygen-oxygen 
and calcium- oxygen interactions are modeled via a Buckingham potential. The 
covalency within the CO3" group is modeled by a Morse potential plus a three-body 
potential and a torsional potential function to restrain the planar CO3" group from 
folding. The interaction between the calcium ions and water is represented by the 
Lennard- Jones (12,6) potential including molecular flexibility 20 using the SPC/Fw 
model 40 . 

Xiao et al. 21 proposed a new forcefield to describe interatomic interactions for the 
bulk aragonite crystal in an aqueous and/or protein environment using the TIP3P 
water model 41 . In comparison to Raiteri's model, Xiao et al. 21 use Lennard- Jones (LJ) 
potentials rather than Buckingham potentials to describe all van der Waals interac- 
tions. They included a LJ potential for the C_C interaction, which was neglected in 
Pavese's work 18 . 

The simulation model of surface follows the approach of Tasker 42 in which the 
crystal consists of two blocks, each comprising two regions that are periodic in two 
dimensions. Region I contains the surface layers and few layers immediately below; 
these atoms are allowed to relax to their mechanical equilibrium. Region II contains 
those atoms further away, which are kept fixed at their bulk equilibrium position and 
represent the rest of the crystal. The size of region 2 is chosen to be thick enough so 
that the atoms represent the effect of the bulk upon the atoms in region 1 . The size of 
region 1 is generally chosen depending upon the convergence of surface energy. This 
quantity is computed using the following expression: 

U s -U b 

where U s is the energy of the surface block of the crystal, U b is the energy of an equal 
number of atoms of the bulk crystal, and A is the surface area. The energies of 
the blocks are essentially the sum of the energies of interaction between all atoms. The 
surface energy is a measure of the thermodynamic stability of a cleavage plane. On the 
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Figure 8 | Equilibrium morphology of hydrous CaC0 3 polymorphs: (a) Ikaite compared to SEM image 39 ; (b) monohydrocalcite. 



other hand there is a widely used alternative criterion named the attachment energy: it 
is the energy associated with the addition of a stoichiometric layer of material onto the 
surface cut. This implies that for any stable material, the attachment energy must be 
exothermic. Therefore, those faces will tend to grow most rapidly. 
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